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Single molecule DNA sequencing via transverse electronic transport using a 
graphene nanopore: A tight-binding approach 
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We present a tight-binding model study of a two-terminal graphene nanopore based device for sequential 
determination of DNA bases. Using Green’s function technique we investigate the changes in electronic 
transport properties of the device due to insertion of a nucleotide into the nanopore created within a zigzag 
graphene nanoribbon. We propose the actual working principle of the device, by setting the bias across 
it to a fixed voltage and translocating a single-stranded DNA through the nanopore and then recording the 
characteristic current signals corresponding to each of the nucleotides at dynamic condition. Our investigation 
shows that this device may turn out to be a very efficient and reliable sequencing device for biomolecules in 
future. 

PACS numbers: 73.23.-b, 73.63.-b, 87.14.gk, 87.15.Pc 


I. INTRODUCTION: 

The idea of nanopore based DNA sequencing has 
emerged from the need to diminish the cost and increase 
the speed of genome sequencing. Traditional marker 
based DNA sequencing needs much more time and also 
quite costly as all the biomolecules targeted for sequen¬ 
tial determination has to undergo fluorescent labelingi^— . 
Dramatic advancements in the DNA sequencing tech¬ 
niques in the last few years have huge impact on genome 
research^^—. As the understanding of genetic components 
of different diseases now become more vivid, it leads to 
the development of new type of medical treatments and 
diagnostics. Genetic research has many other applica¬ 
tions apart from human medicine, such as in security, 
biology and agriculture. But the main problem with the 
current sequencing technology e.g., the well-established 
Sanger method^, is its cost. Sequencing a specific human 
genome to high quality currently costs around $10K de¬ 
pending on the procedures. Whereas the nanopore based 
DNA sequencing offers marker-free single-molecule de¬ 
tection without any sample modification at much lower 
cost ($1000). The initial attempts of nanopore based 
sequencing of single-stranded DNA (ss-DNA) using pro¬ 
tein pore (a-Haemolysin) were made more than a decade 
ago by translocating ss-DNA through the nanopore un¬ 
der an applied bias while reading the changes in the ionic 
current passing through the poreiSrJ^. Detection of the 
signatures of the nucleotides trapped inside a nanopore 
i.e., of a static strand in nanopore is easyi^, whereas the 
problem of detecting the responses of the nucleotides of 
a translocating DNA through nanopore is still a big chal¬ 
lenge, as the speed of translocation is too high for getting 
the necessary current resolution unless they are slowed 
down by chemical modifications with bulky groups^^. 
The early usage of biological nanopore in DNA sequenc- 
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ing is mainly due to the key factors like nanometer dimen¬ 
sion of the pore, detectable ion-current modification, and 
abundance of different processing enzymes and chemicals 
to regulate the speed of DNA translocation in discrete 
step o^^’^^ . However there are some serious disadvantages 
with biological nanopores due to their poor mechanical 
strength!^, and difficulty in integrating them with on- 
chip electronics^^. Solid-state nanopores overcome these 
drawbacks, and also provide some additional advantages 
like possibility of other types of detection, different from 
ionic current^^— . But one big problem remains, the 
synthetic membranes used for detection are more than 
10 nm thick, which means that several nucleotides will 
occupy the nanopore at the same time (distance between 
two consecutive nucleotides is 0.34 nm in double-stranded 
DNA and 0.5 nm in ss-DNA^Si^) and it becomes impos¬ 
sible to perform single-molecule base-specific measure¬ 
ments. Graphene, a single layer allotrope of carbon^, 
provides a new way to develop nanopore based technique 
for the single-molecule sequencing of DNA. Graphene 
is much more superior over the conventional solid-state 
nanopores as the nanopore dimension can be controlled 
to the atomic length scale as well as the thickness of the 
membrane can be of the order of the distance between two 
consecutive nucleotides in a DNA chain. Graphene shows 
better electromechanical properties^. Graphene based 
nanopore also offers several ways of sequential detection 
by means of transverse tunnelling2ii2^, nanoribbon con- 
ductance^^— and also using mechanical deformation of 
DNA^S followed by some experimental realizations^SilL. 

In the present work we report a tight-binding formula¬ 
tion based investigation of a graphene nanopore sequenc¬ 
ing device. The earlier reports available in the litera- 
tur o^^'^^i^^~ — mostly use ab-initio techniques, whereas 
we find that much simpler tight-binding approach pro¬ 
vides a very good physical understanding of the sequenc¬ 
ing device. Using Green’s function approach and Lan- 
dauer formalism we study the local density of states 
(LDOS), conductance and I-V response to detect the se¬ 
quential arrangement of the nucleotides while trapping 
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them inside the nanopore at static condition as well as by 
translocating a ss-DNA through the nanopore. We report 
the performance of the device on the basis of our model 
calculations and show that the characteristic current am¬ 
plitudes for the four different nucleotides are quite dis¬ 
tinct so that they could be detected unambiguously. The 
present study reveals that the graphene nanopore based 
devices could be used for reliable sequential detection of 
the DNA bases following our proposed working principle. 

II. THEORETICAL FORMULATION: 

In order to determine the sequence of the nucleotides 
we use a graphene nanoribbon with zigzag edges with a 
pore at its centre (Fig. [T|). We create the pore in such a 
way that the bipartite two-sublattice nature of graphene 
is maintained. We use semi-infinite zigzag graphene 
nanoribbon (ZGNR) as electrodes^ also, attached to the 
both sides of the device. The active ZGNR system includ¬ 
ing the nanopore (excluding the electrodes) termed as 
device can be represented by the following tight-binding 
Hamiltonian 



FIG. 1. (Color online). Schematic view of a ZGNR nanopore 
with a ss-DNA passing through it. The color dots in ss-DNA 
represent different nucleotides. ZGNR electrodes are shown 
in the figure on the both sides of the device. 


Hzgni - ^ ^ G -f H.C.^ (1) 

i=l 

where e is the site-energy of each carbon atom in ZGNR 
and t is the nearest neighbour hopping amplitude. Ci and 
cj represent electron annihilation and creation operators 
at the ith Wannier state respectively. The Hamiltonian 
of the ZGNR electrodes can also be expressed in the same 
way. The total Hamiltonian of the entire system can be 
written as Htot — Hzgnr T ^electrodes i where 

Htun = T (cqCi -I- cjyCAT+i -I- H.c.^ is the coupling Hamil¬ 
tonian between the active ZGNR device and the elec¬ 
trodes and the coupling between them is defined by the 


hopping integral t. In our calculations we set t=I.O eV 
and e=0 eV. The Hamiltonian of a ss-DNA chain is given 
by 

N 

Hdna= ^ (eiclci + ti,i+ic|ci+i -I- H.c.^ (2) 

i=l 

where is the site energy of the respective nucleotides 
and tiy+i is the hopping amplitude between them. 

The LDOS profiles of the nucleotides trapped at the 
nanopore are defined in terms of the Green’s function as 
p{E,i) = —ilm[Gii(E)] where, G{E) = {E — H + ir])~^ is 
the Green’s function for the entire system including the 
nucleotides with electron energy E as r] ^ 0'^, H = 
Hamiltonian of the entire system, and, Im represents 
the imaginary part. Gii(E) is the diagonal matrix ele¬ 
ment (< i\G{E)\i >) of the Green’s function, \i > be¬ 
ing the Wannier state associated with the trapped nu¬ 
cleotide. To obtain the transmission probability of an 
electron, we first calculate the self-energy of the left 
(right) semi-infinite zgnr electrode following the 

recursive Green’s function technique^i^. The single¬ 
particle retarded Green’s function for the entire system 
at an energy E is given by G’' = [G“]'l' = [E — Hzgnr — 

+ *^]■^ where = ^tun^lfi-^tun is 

the retarded (advanced) self energy of the left (right) 
electrode and is the retarded (advanced) Green’s 

function for the left (right) electrode. The transmission 
probability of an electron with incident energy E is given 
by T{E) = TrpiG’TijG®], where the trace is over the 
entire Hilbert space spanned by the zgnr and the DNA 
nucleotide and ^l{r) = ~ ^L(fl)] i® i'be linewidth 

function. At absolute zero temperature, in the Landauer 
formalism current through the system for an applied bias 

voltage V is given by I{V) = ^ lEF-eV /2 ’ 

where Ep is the Fermi energy of the system. We have 
assumed that there is no charge accumulation within the 
system. 

III. RESULTS: 

For numerical calculations, the site energies of the 
nucleotides are taken as their ionization potentials, ex¬ 
tracted from the ab-initio calculations of Senthilkumar 
et. alM which gives: eG= 8.178, e^= 8.631, ec= 9.722, 
and eT= 9.464, all units are in eV. Then we shift the 
reference point of energy to the average of the ionization 
potentials of the nucleotides which is e=8.995 eV, and 
with respect to this new origin of energy the on-site en¬ 
ergies for bases G, A, C, and T become -0.82 eV, -0.37 
eV, 0.72 eV, and 0.47 eV respectively. The transport be¬ 
haviour of the system would not be effected due to this 
shift of origin of energy. This is equivalent to some ear¬ 
lier reports^^— where the average of ionization potential 
was set as the backbone site-energy. In order to show the 
efficiency and reliability of the proposed device, we first 
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FIG. 2. (Color online). LDOS of the four nucleotides trapped 
at the nanopore. There are four distinct peaks represent¬ 
ing different bases close to their respective characteristic site- 
energies. 


calculate the LDOS. Fig. [5] shows the variation of the 
LDOS at nanopore site where we trap the nucleotides. 
We study this LDOS response of the bases as individ¬ 
ual i.e., we trap the four nucleotides A, T, G, C one 
by one inside the nanopore and study their LDOS pro¬ 
files. The four nucleotides have distinct LDOS patterns 
with peaks around their characteristic site energies. The 
position and height of the peaks are different signifying 
their different electronic structure. This relative differ¬ 
ence of LDOS corresponding to the four nucleotides gives 
a chance to detect them using ARPES technique by trap¬ 
ping them one at a time inside the nanopore. As the be¬ 
haviour of LDOS is mainly controlled by the nucleotides 
and the effects of backbones are not very important^, 
the LDOS measurement technique could be helpful for 
identification purposes of the nucleotides. 

Fig. [3] shows the transmission spectra of the zgnr- 
nanopore device for different cases. The coupling pa¬ 
rameter between the boundary sites of nanopore and the 
nucleotides is taken as r=0.2 eV. The hopping parame¬ 
ter between two identical nucleotides in the DNA chain is 
taken as ti_i+i=0.35 eV and for different bases ti_i+i=0.17 
eV— . For both the intrinsic zgnr and zgnr-nanopore 
the transmission profiles are symmetric with respect to 
zero of energy (Fig. [3a]) and it gets distrayed as we in¬ 
sert a nucleotide (Guanine) inside the pore (Fig. [3b]). 
This is due to the fact that while creating the nanopore 
we have preserved the bipartite two-sublattice symmetry 
of graphene. As soon as a nucleotide is inserted into the 
nanopore this two-sublattice symmetry breaks down and 
the transmission spectrum becomes asymmetric with re¬ 
spect to E=0. Fig. [3b] depicts the relative changes in the 
transmission probabilities for a bare nanopore and a G- 
nanopore (nucleotide Guanine trapped inside the pore). 
There are two distinct peaks in the transmission profile 
for bare nanopore placed symmetrically at -|-ve and -ve 
energies and in presence of Guanine base these peaks 
still remain but get distorted to some extent in asym¬ 
metric manner showing sensitivity of our device to the 
presence of the biomolecule. This feature of bipartite 


two-sublattice symmetry breaking of graphene during ss- 
DNA sequencing has never been reported earlier. 

In Fig. [3c] we present the transmission profiles of the 
nanopore device inserting each of the four nucleotides 
into the nanopore. The difference between the four trans¬ 
mission spectra are much prominent in the -l-ve energy 
region. In Fig. [3d] we zoom in the curves of Fig. [3c] in 
the energy interval 0 to 1.4 eV for better visualization. 
Thymine and Gytosine bases are clearly differentiable in 
the range 0.2 to 0.5 eV with two distinct characteris¬ 
tic peaks in their transmission profiles. In the midway 
between 0.8 to 1.0 eV Adenine (A) and Guanine (G) be¬ 
come differentiable. After that, at higher energy values 
nucleotides become less differentiable as we go far away 
from their characteristic site-energies. 

Fig. [4a] shows the I-V responses of the device for 
bare-nanopore and G-nanopore. The effect of the nu- 
cleobase becomes prominent at considerable bias (inset 
of Fig. [4a]). Fig. [4b] shows the I-V responses of the 
proposed nanopore device in presence of the four nu¬ 
cleotides. All the four bases are distinguishable from 
the I-V response for most of the applied voltage range, 
though curves are found to be maximally separated in 
the range from 0.5 to 1.2 volt (inset of Fig. [4b]). The 
mutual separation of the curves for T, G, G are quite 
large whereas that between A and G is much smaller, 
this might be due to the fact that both A and G belong 
to the same purine group. 

Finally we present the actual working principle of the 
nanopore device. In Fig. [5] we present the current signal 
of the device as we translocate a random ss-DNA through 
nanopore. We take a 30-base long Random ATGC chain, 
translocate it through the nanopore and record the char¬ 
acteristic current signals of each nucleotides. During this 
translocation the bias is kept fixed at 0.54 Volt, the volt¬ 
age that gives maximum separation between the char¬ 
acteristic currents of the four nucleotides. This can be 
found out by studying the current responses of the de¬ 
vice while each and every nucleotide is trapped into the 
nanopore at static condition (Fig. U]). The relative sep¬ 
aration of the signals between pyrimidine and purine 
groups are much larger than that between two bases of 
the same group. The separation between Adenine (A) 
and Thymine (T) is maximum, while that between Ade¬ 
nine (A) and Guanine (G) is quite low as they belong 
to the same purine group. The maximum and minimum 
values of percentage separation {{Imax - Imin)/Imin) of 
the signals are 54.45% and 4.42% respectively, which 
implies that the signals corresponding to the four nu¬ 
cleotides could be differentiated with much more reliabil¬ 
ity than before^i^i^^i^^— . It is worthwhile to mention 
that though we have presented current in arbitrary unit, 
but if we put the numerical values of the various con¬ 
stants e.g., h, h and e, it turns out that the currents are 
of the order of 10/iA, being much higher than earlier re¬ 
ports as well as much greater than the noise level of this 
type of device which being of the order of nA^. These 
high current output promises that the device can perform 
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FIG. 3. (Color online). Transmission probability T(E) with energy for different cases, (a) Comparison between intrinsic zgnr 
and zgnr-nanopore, two-sublattice symmetry being preserved, (b) Curves for a bare nanopore and G-nanopore: breaking of 
two-sublattice symmetry of graphene, (c) Characteristic features of each of the four different nucleotides trapped inside the 
nanopore, (d) Enlarged view of a section of the plot (c) for better visualization. 


sequential detection even under environmental fluctua¬ 
tions without much uncertainty (see Appendix). A very 
recent report on dynamical effects of environment on the 
conductance of graphene nanopore devices^, shows that 
fluctuations of the nucleotides inside the nanopore may 
change the conductance of the devices relying on tunnel¬ 
ing mechanism, though they conclude that these effects 
would not be very important for the devices which rely on 
transverse conductance with larger transmission proba¬ 
bilities. Whereas a study by Krems et in 2009 dealing 
with different types of noises which may occur in actual 
sequencing experiments showed that these environmental 
effects do not strongly influence the current distributions 
and working efficiency of these devices. Though based 
upon these reports we can say that the overall sensitiv¬ 
ity of the device won’t be hampered too much even in 
presence of environmental effects, as the proposed device 
relies on transverse conductance and produces greater 
current output, but there will always be sources of noises 
in actual experimental condition due to environmental 
fluctuations, presence of water and counterions which 
can affect device operation. Being said that, it is also 
important to note that the results given in this work is 
open to improvements by functionalization of the edge 
atoms of the nanopore which can significantly enhance 
nucleobase-pore interaction, thus reducing the structural 
noise by enhancing the graphene-nucleobase electronic 


coupling^SiS. Different types of groups can be used for 
functionalization (e.g., hydroxide^, amine or nitrogen^) 
to provide custom made solution to overcome noise in 
electrical DNA sequencing techniques. 


IV. CONCLUSION: 

We propose a reliable DNA sequencing technique us¬ 
ing graphene-based nanopore device on the basis of our 
model calculations. Most of the previous proposals on 
graphene nanopore based detection techniques rely on 
tunnelling currents that are mediated via energy levels in¬ 
troduced by the nucleotides^Mii^^iSI. On the other hand 
the current device relies on comparatively large trans¬ 
verse conduction current which gets modulated in pres¬ 
ence of the DNA bases, and hence gives much high resolv¬ 
ing power. Presence of a nucleotide inside the nanopore 
modifies the LDOS at and around their characteristic 
site-energy and affects the transmission profiles only in 
the neighbourhood of this energy. The resulting trans¬ 
verse current is also quite sensitive to DNA bases, and 
hence it renders graphene-nanopore system as a very ef¬ 
ficient sequencing device for biomolecular detection. We 
also propose an unique way of operation of the device, by 
fixing the bias voltage across the device and translocat¬ 
ing a ss-DNA through the nanopore, while recording the 
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FIG. 4. (Color online). Current - Voltage response of the ac¬ 
tive nanopore device for different cases, (a) Curves for a bare 
nanopore and G-nanopore. (b) Characteristic current out¬ 
puts for each of the four different nucleotides trapped inside 
the nanopore. The insets show the enlarged view at important 
voltage ranges. We set the Fermi energy at Ep—Q eV 



Base index 

FIG. 5. Stop and go translocation of a Random ATGC ss- 
DNA chain through the nanopore, while bias across the device 
is fixed to a specific value (0.54 Volt) which gives maximum 
separation in current responses for different nucleotides. We 
record the characteristic current output for the nucleotides 
as they translocate through the nanopore. Respective nu¬ 
cleotides are represented within the figure by their usual sym¬ 
bols A, T, G and C adjacent to their characteristic current 
outputs. 


distinct current responses corresponding to the different 
nucleotides. With the two-sublattice symmetry breaking 
of graphene due to insertion DNA nucleotides, large cur¬ 
rent output and appreciable amount of accuracy, we have 
been able to put forward an unique operating principle 


for sequential detection of DNA bases that also gives un¬ 
matched accuracy even under environmnetal fluctuations 
(see Appendix) compared to previous results2ii2ii2&iiiivi&. 
In summary, our proposed technique and model calcula¬ 
tions provide a comprehensive details of the effects of 
the four nucleotides. Hence we think, following this work 
high speed unambiguous sequential detection of the DNA 
bases can be achieved at a much lower cost. 


V. APPENDIX: 

In this section we provide some additional calculations 
on sequencing of ss-DNA. In Fig. |6]we present the vari¬ 
ation of current responses of the device for a translocat¬ 
ing ss-DNA chain as we incorporate the fluctuations of 
the nucleotides inside the nanopore during translocation. 
Due to these fluctuations DNA nucleotides will undergo 
several random orientations inside nanopore and it will 
result in changes of nucleotide-to-nanopore coupling^^. 
We simulate these fluctuations in our model by consid¬ 
ering the coupling between edge atoms of nanopore and 
the nucleotides (t) to be randomly distributed within the 
range 0.01 to 0.2 eV. It is clear form the figure that due 
to these fluctuations the magnitude of the responses get 



FIG. 6. Stop and go translocation of a 30-base long ran¬ 
dom ATGC ss-DNA chain through the nanopore, incorporat¬ 
ing the fluctuations of the nucleotides within the nanopore 
during translocation while bias across the device is fixed at 
0.54 Volt. We record the characteristic current output for the 
nucleotides as they translocate through the nanopore. The 
sequence of the DNA chain is GTCCAGAGAAGCTACACG- 
GATTTTATCCGT. 

decreased from their characteristic values (Fig.5) and 
there is also variation in the characteristic current am¬ 
plitudes. As an example current output for G varies 
from 0.1078 to 0.1097 (arb. unit). For T current am¬ 
plitude varies from 0.1260 to 0.1337 in arb. unit which 
has decreased from its characteristic value of 0.1590 (arb. 
unit). Similarly for C current output varies within the 
range 0.1200 to 0.1219, decreased from 0.1288 and for 
A current varies from 0.1023 to 0.1059, changed from 
characteristic value of 0.1029, all units are in arb. unit. 
From these specific values we can see that variations are 
not very large and they are not at all overlapping with 
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each other, which shows that even in presence of this 
kind of environmental effects our device can work with¬ 
out much ambiguity. Even the drops in their character¬ 
istic current outputs are also small, maximum drop has 
been observed for T and minimum for A (in some cases 
current increases). The lowest value of current output 
achieved under these fluctuations is 0.1023 in arb. unit 
(that is actually 0.1023x10=1.023 /xA) which also shows 
that the device can produce higher current output even 
in presence of this type of fluctuations. 
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